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Abstract 

Based on dynamical renormalization group (RG) calculations to the one- 
loop order, the surface growth described by a nonlinear stochastic conserved 
growth equation, ^ = ibz^2V^/i + AV-(V/i)'^ -\- rj {v2 > 0), is studied ana- 
^ . lytically. The universality class of the growth described by the above equa- 



tion with +1/2 (diffusion) is shown to be the same as that described by the 



o . 

CJ ■ Edwards- Wilkinson (EW) equation (i.e. +V2 and A = 0). In contrast our RG 



• ^ . recursion relations manifest that the growth described by the above equation 

>< 

^ . with —1/2 (anti-diffusion) is an unstable growth and do not reproduce the re- 

cent results from a numerical simulation by J. M. Kim [Phys. Rev. E 52, 
6267 (1995)]. 
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Recently, there have been considerable interests in the various surface growth models |l| 
to understand roughenings in the growing surfaces theoretically. Since the surface structures 
of many growth process are self-affine, most efforts have concentrated on the surface width 
W, which is defined by the root mean square fluctuation of the surface height. In a finite 
system of lateral size L, the width W starting from a fiat substrate scales as 

W{t) ~ L^it/L^ 

r-t", t < (1) 

~ L", t > L' 

where the scaling function /(x) is for x <^ 1 and constant for x 3> 1. The exponents (3 
and z are connected by the relation z = a/13. Among the growth models which satify Eq. 
(1), the class of models known as "conserved" models [0-^, which conserve the total number 
of particles after being deposited, has been extensively studied as a possible description for 
the real molecular beam epitaxy (MBE) growth ||T0|. In these conserved models the height 
h describing the local position of the surface of growing materials is known to obey 

^^ = -V-J(x,t) + ry(x,t) , (2) 

where J is the surface current and r7(x, t) is an uncorrelated random noise 

(r/(x, t)r7(x', t')) = 2D5\^ - ^)5{t - t') . (3) 

There have been several studies on the surface growth described by the following current 
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J(x,t) = -i^2V/i + z/4V^/i- A(V/i)^ (4) 
and thus by the corresponding continuum equation 

dh{^ t) ^ ^^^2^ _ ^^^4^ ^ AV- (V/i)^ + r/(x, t) . (5) 

For 1^2 > 0, the critical property of the surface growth described by Eq. is solely 
dependent on the first term 1'2'V'^h in the right hand side, which is the most relevant term 



in a renormalization group (RG) sense and it belongs to the Edwards and Wilkinson(EW) 
universality class with a = {2 — d)/2 and z = 2 for any substrate dimension d ||Tl],|12|. 



When = A = 0, Eq. (|^) becomes the MuUins-Herring equation [|I^, which can be solved 
exactly to give a = {4 — d)/2 and z = 4 0. When 1^2 = 0, z/4 > and A > 0, the recent 
studies |n],|T^ have shown that the most relevant nonlinear A term generates an effective 



h'2^-^V'^h term with 1^2^^ > ^he corresponding growth belongs to the EW universality 
class, even though both the scaling argument and the dimensional analysis 0] suggest the 
exponent values /3 = (4 — d)/2{A + d) and a = (4 — (i)/4. 

More recently the growth described by Eq. (5) with U2 < 0, > and A > has been 
studied numerically [jl3[. From this numerical results it has been argued that the critical 
property of such growth belongs to the EW universality class even for 1/2 < 0, on the ground 
that the nonlinear A term could supress any negative i>2, so that 1^2 could become effectively 
positive, i.e. 1^2^^ > IHI- If this argument is right, then the A term could make the unstable 



growth with -|z/2|V^/i term the stable EW growth with 1/2^^ > 0. However if one analyzes 
Eq. from a point of view of the scaling and dimensional analyses, the most relevant term 
is 1'2'V'^h regardless of the sign of z/2. So it is somewhat strange that the less relevant A term 
renormalizes the most relevant z/2 term so effectively as to make the change of sign of z/2. It 
is thus our motivation to check whether this strange renormalization is physically plausible 



or not by use of dynamical renormalization group calculations |19|J2^ and other analytical 
methods. 

To achieve this goal, we first discuss the dynamical renomalization group study of Eq. 
(^) with z/2 7^ 0, z/4 = and A > 0. In the large distance and long-time hydrodynamic limits, 
Eq. for z/4 = in Fourier space can be written as 

h{\i,uj) = G(k,Lj)ri(k,Lj) 

= Go{k,uj)r]{k,u) + XGo{k,uj) 
dQid'^qidn2d'^q2 



(27r)2rf+2 

x(qi ■ q2)[k • (k-qi-q2)]/i(qi,r]i) 



x/i(q2, ri2)/i(k-qi-q2,cj-f2i-f22) 



(6) 



where Go(k, cu) is the bare propagator defined by the expression 
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Go(k,cu) 



(7) 



Using {rjCk, uj)ri{k' , uj')) = 2D6'^(k + k')5(co' + uj') and performing the internal frequency 

integrals of Eq. (|]) to the one-loop order, one obtains 

XD d + 2 



Gik,iu) = GoiKcu) K,—-k'Gt,{k,uj)J:iq), 

i>2 d 



(8) 



in the limits u ^ and A; 0. Here = dqq'^ ^, A is the momentum cutoff (A = 1), 
and Kd = Sd/{2nY with Sd = 27^'^/'^ /T{d/ 2). Since has no infrared divergence for any 
dimension, we can expect that the nonlinear A term is irrelevant in RG sense and that the 
A will not renormalize negative z/2 term to be positive. 

For the confirmation's sake, we want to calculate the dynamical RG recursion relations 
for V2 and A. Integrating out the fast modes in the momentum shell G [Ae^^,A] and 
restoring the slow modes k"^ G [0, Ae^^] with = 1 — (5^ + - , we find an effective surface 
tension z/^ for the long wavelength modes, 

XDd + 2 



In an analogous way we get 



V2 



A< = A 



1 + 5iKd 



vl d 



1 - 6iKda{d)^ 

^2 



(9) 



(10) 



where a{d) = 9 for d = 1, and a{d) 



for d = 2. Upon requiring that the equation 



stays invariant under the scale(e ) transformations k e~ k, t e t and h e , the 
parameters should transform as follows: z/2 6^"^ z/2, D ]jz~2a~djj _^ ^z+2a-4^_ 

Combining these scale transformations and Eqs. (P) and (|10]), we get the following RG 
recursion relations: 



dV2 

Hi 

dD 

~di 

dX 



1^2 



z~2 + gKd 



d + 2 
d 



D[z-2a-d] 
^ = X[z + 2a-4:- ga{d)Kd] , 



(11) 
(12) 
(13) 



where the couphng constant g is defined as 



(14) 



The RG recursion relation of the effective couphng constant g is thus 

'^^=g[-d-gA{d)K,] , (15) 

where ^(1) = 15 and A{2) = y for d = 1 and 2 and we have also confirmed A{d) > for 
any d. As one can expect from Fig. 1, the RG fiow of g has only one stable fixed point which 
is a Gaussian fixed point, i.e. g* = and does not have any nontrivial fixed points. The 
RG fiow of g for z/2 < is the same as those for z/2 > 0. If the results from the numerical 



calculations in Ref. [jT3[ would be explained by the dynamical RG, then the RG fiow of g for 
z/2 < should be different from those for z/2 > 0. As shown in Fig. 2, the RG flows starting 
at (z/2 > 0, A > 0) will eventually arrive at {1^2^^ > 0, A — 0). This means that the growth 
described by Eq. (5) with z/4 = 0, z/2 > and A > in the large distance scale belongs to 
the EW universality class as we have expected. The RG flows starting at (z/2 < 0, A > 0), on 
the while, will eventually arrive at (z/2^^ < 0,A ^ 0). This means that the growth described 
by Eq. (5) with z/4 = 0, z/2 < and A > in the large distance scale has nearly the same 
critical property as the unstable linear growth described by 

dh{^,t) 



dt 



-\u2\V'h + r] . (16) 



Based on our dynamical RG calculations we conclude that the growth by Eq. (5) with 
z/2 < 0, A > and z/4 = should be an unstable growth and we cannot reproduce the recent 
results by a numerical simulation ||T3| . 

The RG recursion relations (11), (12) and (13) are based on the one-loop order pertur- 
bative calculations and are intrinsically exact only for A^loro = ^<^l. The main 
result of the dynamical RG calculations physically indicates that the growth described by 
Eq. (5) with z/2 < 0, z/4 = and the small A(> 0) does not belong to the same universality 
class as that of the EW growth. To understand physically the unstable growth by Eq. (5) 



with 1/2 < 0, A > and 1/4 = for the finite A including the case for A ^ 1, let's discuss the 
growth from an Hamiltonian-based argument or from an equlilibrium physics. If one believe 
that the dynamical equation (|^) with 1/4 = can be derived by an equilibrium Hamilto- 
nian Ti via the Langevin equation ^ = —^ + V and the Langevin equation will reach a 
steady state or an equilibrium where the states are controlled by the Boltzman distribution 
PiH) oc exp{—(3T-C), H should be 

n = jd'^xqivhf + \{vhY] . (17) 

If one puts Vh to be equal to a field (i.e., Vh = (f){x) ) and use a mean- field theoretic 
argument for Hamiltonian (|l^, the corresponding Landau- Ginzburg function C becomes 

C = "-1^^ + . (18) 

When z/2 > and A > 0, the mean field theory with |4 = gives (</)) = (V/i) = as in a 



disordered phase of the ordinary magnetic phase transitions pT|. The fiuctuations around 
the mean-field (0) = (V/i) = should be a Gaussian-type and it seems quite plausible to 
believe that the surface roughenings with z/2 > and A > belong to the EW universality 
class. However when z/2 < and A > 0, the mean field thory predicts 

(0') = ((V/.)^)e, = -y > (19) 
as in an ordered phase of the ordinary magnetic phase transitions. In the surface growth phe- 



nomena, Eq. ([19|) means that the local slope satisfies (|V/i|)eg — y |z^2|/A in an equilibrium 
state, i.e. in a saturation state. Dynamically, in the growing process when (|V/;-|) < y^|z/2|/A, 
the anti-diffusion term — Z/2V/1 of the current (see Eq. (^)) dominates. Dynamical RG re- 
cursion relations (11), (12) and (13) which are intrinsically exact for A^lor5f = ^^l 
correspond to this growing stage of the local slopes. For (7^1 where RG calculations are 



correct, it takes very long time for the local slopes to reach the equilibrium value a/|z/2|/A and 



this dynamical process should be described by Eq. (16), because {\Wh\)eq — y|z^2|/A ^ 1. 
We have confirmed this kind of growth by using several simulations for the growthes with dif- 



ferent gf's [^. In contrast for finite g^s the time interval in which the anti-diffusion dominates 



become rather finite and sooner or later the growth reaches the saturation state or the equi- 



librium state where the effect of the anti-diffusion term is balanced to give (|V/i|) ~ y |f2|/A. 
After that local slopes flucuates around the value in Eq. (19). We have also confirmed that 



this kind of the unstable growth occurs when g is finite by simulations |2^. It is thus not 
physically sound that the growth with z/2 < and A > belongs to the same universality 
class as that of the EW-like growth with with i>2^^ > as in Ref |T^. Instead we 

do expect an unstable growth. 



We now want to discuss on the point why the numerical simulations in Ref. |T3[ couldn't 



see such an unstable growth and did see only the stable growth which belongs to the EW 



universality class. To do a numerical simulation in Ref. ||T3l, the local currents in Eq. @ in 
d = 1 has been set to be equal to 

j(k,i) = U2[h{k) - h{i)] - Ui[h{k + 1) + h{k - 1) 
-2h{k) - h{i + 1) - h{i - 1) + 2h{i)] 

+X[hik) - hit)f , (20) 

where k is either i + 1 and i — 1 and a particle is added on site k if j{k) is negative. If both 
j{i + and j{i — are positive, a particle is added on site i. In case both j{i + 
and j{i — are negative, a particle is added on either site randomly. When 1/4 = 0, Eq. 
(pO|) can be written as 

J{k,^) = U2[h{k) - + -im - h(z)]'} . (21) 

In the initial growing process with z/2 < when {[h{k) — /i(i)]^) is quite small and { } term 
in Eq. ( pT]) is positive, the anti-diffusion should dominate the growth. As the growth is 
developed, {[h{k) — h^i)]"^) increases to —h'2/X rapidly and after that { } term in Eq. (|2lD 
fluctuates around zero. We have conflrmed that at this stage the local slope {\h{k) — h{i)\) 



is nearly equal to y— 2/2/^ In contrast the numerical simulations in Ref. ||13[ have 

been done only for the cases —U2l\ = 1,3 and the corresponding saturated local slopes 
\h{k) — h{i) \ are only 1 and v^- Since h{k) must be integer numbers 0,1,2,... in a simulation. 



\h{k) — h{i)\ = 0, 1 or 2 could be easily generated even by the noise effect only, so the { 
} term in Eq. ( pT]) being negative for —1'2/X = 1,3. This fact explains why the numerical 
simulation in Ref. [0 had seen only the EW-type growthes instead of the unstable growth. 
But for the considerably large —1'2/X, the anti-diffusion term dominates from the beginning 
and this effect stays longer and longer until being cancelled out by the nonlinear A term. 
Hence hardly realizes the EW-like behavior. 
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FIGURES 

FIG. 1. The RG flows of the couphng constant g. There is only one fixed point g* = which 

is attractive fixed point. The RG flows of g for f2 < are the same as those for 1/2 > 0. (See Eq. 
(15)) 

FIG. 2. Schematic RG flows in (A, f2)-plane. for d = 1 and 2, respectively. The RG flows 
starting at (1/2 > 0, A > 0) will eventually arrive at {i'2'^^ > 0,\ —>■ 0). In contrast the RG flows 
starting at (^2 < 0, A > 0) will eventually arrive at {1^2^^ < 0, A — > 0). 
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